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1 Introduction 



A fundamental test of QCD is obtained by comparing the value of the quantum chromody- 
namics (QCD) coupling constant a$ extracted from experimental data at different energy 
scales with theoretical predictions given by the renormalization group method. The corre- 
sponding evolution law of as{Q 2 ) is now experimentally studied down to low momentum 
transfers, Q 2 ~ lGeV 2 . The principal experimental information concerning as at low Q 2 
comes from measurements of the inclusive hadronic r-decay rate and from deep-inelastic 
scattering sum rules. As is well-known, theoretical QCD predictions for these processes 
are primarily based on perturbation theory (PT) improved by the renormalization group. 
However, the perturbative expansion is ill-defined at low energies and the conventional PT 
method of deriving the QCD running coupling constant leads to unphysical singularities of 
&s{Q 2 ), such as a ghost pole, which are in conflict with the fundamental principle of causal- 
ity. Higher-order PT corrections taken in the asymptotic form cannot resolve this problem 
and just add unphysical cuts. Moreover, at low values of the momentum transfer there 
is a strong dependence on the choice of renormalization scheme, which leads to essential 
ambiguities in the description of the physical quantity under consideration. 

In IJ, the conventional method of the renormalization group improvement of the PT 
approximation has been modified by requiring Kallen-Lehmann analyticity, which reflects 
the principle of causality. As a result, the QCD running coupling constant has no unphysical 
singularities and at low energy scales behaves significantly differently than the conventional 
perturbative one. The method proposed in [|TJ, called analytic perturbation theory (APT), 
gives a well-defined procedure of analytic continuation of the running coupling constant from 
the Euclidean (spacelike) into the Minkowskian (timelike) region [0, [|. This fact allows one 
to describe the inclusive decay of the r lepton in a self-consistent way Q, essential became 
the correct analytic properties are important in order to rewrite an initial integral expression 
for the i? T -ratio in the form of a contour integral representation. 

The renormalization scheme (RS) ambiguity which appears due to the truncation of the 
perturbative expansion represents a serious difficulty in carrying out perturbative calcula- 
tions. There is no definite solution to this problem apart from computing indefinitely many 
terms in the perturbative expansions. To somehow avoid this difficulty, various optimization 
procedures have been applied, for example, the principle of minimal sensitivity || or the 
effective charge approach 0. The RS dependence can be also reduced by using the Pade 
summation method 0j. As it has been argued in f| ||, besides considering PT predictions 
in some preferred scheme, one should also study the stability of those predictions by varying 
parameters that define the RS over some acceptable interval. In the framework of the ana- 
lytic approach the problem of the RS dependence has been studied for the process of e + e~ 
annihilation into hadrons and the process of inclusive r-decay in (1^, [11]], where it has been 
demonstrated that the APT approach can reduce the RS dependence drastically. 



Apart from the r decay, there is another important observable, the Bjorken sum rule [12 
which allows one to extract a$ at low Q 2 . In this paper we analyze the Bjorken sum rule, 
which is the integral of the difference between gi for the proton and neutron, in the framework 
of the APT method. We compare our result with standard PT predictions (see, e.g. |jl3| ) 
and investigate the RS dependence of the various theoretical predictions. 
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2 The Bjorken sum rule within the analytic approach 



The polarized Bjorken sum rule refers to the integral over all x at fixed Q 2 of the difference 
between polarized structure functions of the proton g\ and the neutron gf, 



rr n (Q 5 



gl(x,Q 2 )-gUx,Q 2 



dx . 



(1) 



The Bjorken integral (H) can be written in terms of the QCD correction A B j 



rr n (Q 2 
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1 - A Bj (Q 2 
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The value of the nucleon beta decay constant is taken to be \gA/gv\ — 1-2601 ± 0.0025 |14 



The perturbative QCD correction to the Bjorken sum rule in the three-loop approxima- 
tion with the use of the MS renormalization scheme and in the massless quark limit has the 
form 
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where for three active quarks the coefficients are d^ s = 3.5833 and d^ s = 20.2153 \TE . The 
perturbative running coupling constant «pt(Q 2 ) is obtained by integration of the renormal- 
ization group equation with the three-loop /3-function. 

As it has been demonstrated in |16| by using the Deser-Gilbert-Sudershan representation 
for the virtual forward Compton amplitude, the moments of the structure functions are 
analytic functions of Q 2 in the complex Q 2 -plane with a cut along the negative part of the 
real axis. It is clear that the perturbative representation (|J) violates these analytic properties 
due to the unphysical singularities of the PT running coupling constant for Q 2 > . To 
avoid this problem we apply the APT method, which gives the possibility of combining the 
renormalization group resummation with correct analytical properties of the QCD correction 
to the Bjorken sum rule. 

Let us write down the QCD correction in the form of a spectral representation 



A Bj (Q 2 
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o + Q 2 
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where we have introduced the spectral function, which is defined as the discontinuity of 
A Bj (Q 2 ): Q(a) = Disc{A Bj (-a-ie)}/2i. 

By calculating the spectral function g(a) perturbatively we get a expression for A B j(<5 2 ) 
which has the correct analytic properties and therefore has no unphysical singularities. For 
instance, the APT running coupling constant in the one-loop approximation has two terms: 



«apt(<5 2 ) 
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Q 2 /A 2 



(5) 



Obviously, the first term in Eq. has the standard PT form, but the second term (which 
appears automatically from the spectral representation and restores the correct analytic 
properties of the running coupling constant) has an essentially nonperturbative nature. If 
we rewrite the second term in Eq. (Rf) in terms of the PT coupling constant, we obtain an 
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expression which has an essential singularity like exp (— 47r//5 apT) as apx — > . Therefore, 
the second term in Eq. (^) does not contribute to the conventional perturbative expansion. 
It has been argued in |] that a similar situation holds also for the running coupling constant 
in higher order approximations. The asymptotic PT expression for the running coupling 
constant is an expansion in the small parameter l/ln(Q 2 /A 2 ). This approximation violates 
the Q 2 -analyticity of the running coupling constant and does not allow one to describe low 
energy scales - the perturbative series diverges in the infrared region. The APT method 
removes this difficulty and leads to a quite stable result for the entire interval of momentum. 
The difference between the shapes of the PT and APT running coupling constants becomes 
significant at low Q 2 -scales. 

It is convenient to write the three-loop APT approximation to ABj(Q 2 ) as follows 

Kf T (Q 2 ) = $apt(Q 2 ) + d^UQ 2 ) + d 2 5$ T (Q 2 ) , (6) 

where the coefficients d\ and d 2 are the same as in Eq. (^) and the functions <5apt(Q 2 ) are 
derived from the spectral representation and correspond to the discontinuity of the fc-th 
power of the PT running coupling constant 

«t«? 2 ) = 1^ jfgi I- Kt(-o - *)} ■ (7) 

The function S^p T (Q 2 ) defines the APT running coupling constant, oaptIQ 2 ) = ^apt(Q 2 ) 5 
which in the one-loop order is given by Eq. fl5|). 

In the case of PT the QCD correction to the Bjorken sum rule is represented in the form 
of a power series in a PT [see Eq. @], but in the case of APT the same QCD correction is 
not a polynomial in the APT running coupling constant. As follows from Eq. (P), the first 
term of the expansion is a apt /n, but the following terms are not representable as powers 
of a apt • 

Table 1: Terms of PT and APT series for the QCD correction to the Bjorken sum rule. 



rr n (3GeV") = 0.160 




IsMerm 2nd-term 


3rd-term 


APT _ 


0.131 + 0.062 + 


0.045 


A APT 
A B.j 


0.190 + 0.045 + 


0.003 



To illustrate the difference between the convergence properties of the PT expansion @ 
and the APT series (§) we use, as an example, a typical value of Ti~ n = 0.160 ± 0.014 
at Q 2 = 3 GeV 2 taken from 0, fixing in such a way the value of the QCD correction in 
Eq. (Q). In Table [I] we present numerical results for contributions to Aej in different orders 
for both the PT and the APT methods. From this table one can see that the higher order 
corrections to the Bjorken sum rule play a different role in the PT and APT approaches. 
The convergence of the APT series seems to be much more well behaved than is that of the 
PT expansion at such small Q ~ 1.73 GeV. 

Next in Fig. [I], we show T^ -11 as a function of the QCD running coupling constant as 
in the PT and APT approaches. As outlined above in the PT case, the function T^ -11 is an 
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explicit function of the PT running coupling constant and in the one-loop approximation is 
represented by a straight line in Fig. [j], as a parabola in the two-loop case, and as a cubic 
curve in the three-loop one. At sufficiently large values of as ~ 0.5, the difference between 
the 1-, 2-, and 3-loop PT predictions becomes very large. In the case of APT, the function 
r^ 11 only in one-loop approximation is an explicit function of the APT running coupling 
constant. Of course, the coincidence of the one-loop PT and APT curves in Fig. [I] does not 
mean that the PT and APT approaches are physically identical, this is simply a matter of 
the linear form of the one-loop approximation - because the behavior of the PT and APT 
running coupling constants are rather different [see Eq. (§)] . The contribution of the higher 
loop corrections in the APT case is not so large as in the PT one and the corresponding 
curves in Fig. [I] are quite close to the linear function. Fig. [I] demonstrates that the APT 
result is more stable with respect to higher loop contributions. 



Table 2: The QCD parameters extracted from different experimental inputs. 



Input [Ref.] 


Experiment 


QCD parameters 




Ql = 10 GeV 2 


rrW Abj(Qo) 


A PT a PT (Ql) A A pt 


ttAPT(Qo) 


(a) [13 


0.183 ±0.034 0.129^0.162 


467 MeV 0.275 741 Me V 


0.301 


(b) 


0.195 ±0.029 0.072^0.138 


138 MeV 0.177 149 MeV 


0.179 



To compare the Q 2 -evolution of the Bjorken sum rule in the APT and PT approaches, 
we use as input new experimental values of Ti~ n (Ql = 10 GeV 2 ) given by the SMC Col- 
laboration in Refs. W7\ and [|T1|. These values are presented in Table |^ as input (a) and 



(b). The experimental errors of the present data are too large to determine the QCD cor- 
rection. However, to illustrate the evolution, we fix the values of the parameters A or the 
corresponding values of as in the PT and APT cases from the central values obtained from 
the data (a) and (b). By using these normalization points, we plot, in Fig. |2|, the Q 2 de- 
pendence of the QCD correction Apj. Besides the normalization points, for illustration, we 
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Figure 2: The 3-loop QCD correction to Bjorken sum rule in the APT and PT approaches 
at different experimental normalization points. 



also represent in Fig. |] the recent data of the E154 Collaboration [19| obtained for small 
Q 2 . Normalizing at the point (b) we get small values of A which are close to each other 
in the PT and APT cases (see Table 0); therefore, the corresponding curves in Fig. |2| are 
also close to each other. However, the values of the running coupling constant extracted in 
such a way are too small to have good agreement with other experimental data, for example, 
with the data of the E154 Collaboration plotted in Fig. ^] and with the value of as extracted 
from the semileptonic r decay ||14|| . If the more realistic normalization at the point (a) is 
used, the difference between the PT and APT predictions becomes large. Now, we have the 
value Aapt = 741 MeV [see Table input (a)] that is consistent with value extracted from 
^ pT = 871 ± 155 MeV. The discrepancy of the PT value A PT = 467 MeV and 



r decay fTJ, A 

that obtained from the r decay, A PT = 385 ± 27 MeV, appears to be large. 

The dotted curve in Fig. |] demonstrates that the inclusion of the higher-twist (HT) 
corrections does not change the overall picture of the Q 2 -dependence of the PT prediction in 
the interval under consideration. However, the higher-twist corrections do reduce the value 
of A. By using the higher-twist coefficient cht = —0.03 GeV 2 (see discussion in |I3|), we 
obtain the value of Apt = 387 MeV which agrees well with the above value of Ap T . Inclusion 
of such higher-twist corrections is, however, not necessary in the APT approach to achieve 
agreement between the QCD scale parameter A extracted from the Bjorken sum rule and 
r decay. Unfortunately, the experimental errors on the Bjorken sum rule are too large to 
reach a definite conclusion. 
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3 Renormalization scheme dependence 



A truncation of the perturbative expansion leads to some uncertainties in theoretical predic- 
tions for a physical quantity arising from the RS dependence of the partial sum of the series. 
At low momentum scales these uncertainties may become to be very large (see, for example, 
an analysis in |20|| ). Thus, it is not enough to obtain a result in some preferred scheme, but 



rather it is important to investigate its stability by varying the parameters that define the 
RS over some acceptable domain. 

Consider the RS dependence of our results. The coefficients d\ and d 2 in Eq. (|3|) are RS 
dependent. In the three- loop /3-function 

d a 

(3(a) = /U 2 7r^7 = ~Po a 2 (1 + cia + c 2 a 2 ) , (8) 
o fl z 

the coefficient c 2 also depends on RS. 

In the framework of the conventional approach, there is no resolution of this problem 
of the RS dependence apart from calculating indefinitely many further terms in the PT 
expansion, and there is no fundamental principle upon which one can choose one or another 
preferable RS. However, it is possible to define a class of 'natural' RS's by using the so-called 
cancellation index criterion |J. According to this criterion a class of 'well-behaved' RS's are 
defined a such a way that the degree of cancellation between the different terms in the second 
RS-invariant 

p 2 = d 2 + c 2 - d\ - di c\ (9) 

is not to be very large. The degree of these cancellations can be measured by the cancellation 
index || 

C = -i-r (\d 2 \ + \c 2 \ +d\ + \dt\cx) . (10) 

\P2\ V ' 

By taking some maximal value of the cancellation index C max one should investigate the 
stability of predictions for the RS's with C < C max . In the case of the MS-scheme, for the 
Bjorken sum rule, the value of the cancellation index is C^g = 8 . We will consider this value 
boundary for the class of 'natural' schemes. 

Three RS dependent parameters d\ , d 2 and c 2 are connected by the second RS-invariant @ 
and, therefore, any RS is defined by the pair of numbers (di, c 2 ). In Fig. |3] we plot the PT pre- 
dictions represented by dashed lines for three schemes: MS; A with parameters d± = —4.3 , 
c 2 = , and with the value of cancellation index C\ = 8 ; and B with parameters df = , 
cf = 14.5 , and Cb = 4.3 . As a example of an optimal scheme, we consider the effective 
charge (ECH) approach in which df CH = d 2 CH = , cf CH = p 2 , and C E ch = 1 • The dotted 
curve presents the ECH result. For the sake of illustration we also show the experimental 
data taken from |T7], H5, |2T|.Q To normalize all curves, we used the experimental value for 



Abj at Ql = 10 GeV 2 given in Table || as input (a).0 The description of experimental data 
within the MS scheme seems to be quite good; however, as has been mentioned above, there 
is no reason why the MS scheme is preferable over, say, the A-scheme. The corresponding 



1 If one returns to Fig. [l] and compares the properties of convergence of the PT and APT series for the 
data pH , where T^~ n (Q 2 ~ 1 GeV 2 ) ~ 0.1 , we can see that the PT series is not stable with respect to higher 
loop corrections. 

2 To this end, we used C < 8 for the B-scheme in order to satisfy the normalization condition. 
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Figure 3: Renormalization scheme dependence of predictions for Abj vs. Q 2 for the APT 
and PT expansions. The solid curves, which are very close to each other, correspond to the 
APT approach in the ~MS, A, B and ECH schemes. The PT evolution in ~MS, A and B 
schemes (dash), and in ECH (dot) are shown, as are the PA results in MS, A, B schemes 
(dash-dot). The SMC data JT7|/ is denoted by a square, the triangle is the E154 data jTQ], 
and circles are E143 data 
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series have the forms A^ s = x + 3.583a; 2 + 20.215a; 3 and — x — 4.36a; 2 + 16.834a; 3 , where 
x = a s (Q 2 )/7r . 

Thus the uncertainties coming from the RS dependence of perturbative calculations are 
rather large. At the same time, the APT predictions (solid curves) practically coincide with 
each other and, therefore, are RS independent. This stability reflects the existence of a 
universal limiting value U of the APT running coupling constant and the small values of 
higher loop corrections for the entire interval of momentum. We evaluated also Ab]{Q 2 ) 
using the Pade approximant (PA) [0/2] of the PT series (see f7|). The results are shown 
as three dash-dot curves corresponding to the MS, A and B schemes. The PA improves 
the stability properties, but the sensitivity to the choice of RS becomes very large for small 
momentum, Q 2 < 5 GeV 2 . 

Thus the conventional PT prediction at small momentum transfers has a very large RS 
ambiguity. The APT approach reduces the RS dependence drastically. At the three loop 
level the APT result is, practically, RS independent. 



4 Summary and conclusion 

We have considered the Bjorken sum rule by using the APT approach and have demon- 
strated that the convergence properties of the APT are much better than are those of the 
PT expansion. The APT results have extraordinary stability with respect to higher loop 
corrections and also to the choice of the RS. The analysis performed shows a quite different 
Q 2 evolution of the Bjorken sum rule in the PT and APT descriptions. At low Q 2 of order 
a few GeV 2 , the conventional PT approach leads to a very rapidly changing function for 
the Bjorken integral in many RS's. At the same time, the three-loop APT prediction is 
practically RS independent and the Q 2 evolution is described by a slowly changing function. 
Unfortunately, the present experimental data for the Bjorken sum rule have large errors, 
which does not allow us to discriminate between the approaches experimentally. Neverthe- 
less note that experimental data of the E143 Collaboration, which have just appeared |2"2" 



show, practically, the same value of Ti~ n at Q 2 = 2, 3 and 5 GeV 2 , which agrees well with 
the slow APT Q 2 evolution. 

The APT approach incorporates "perturbative" power corrections to secure the required 
analytic properties of the running coupling constant. These corrections come from the per- 
turbative short distance analysis and their appearance is not inconsistent with the operator 
product expansion |p3| . In this note we have concentrated on the analytically improved 
perturbative contribution to the Bjorken sum rule and did not consider higher-twist effects. 
The higher-twist terms can, in principle, be potentially important; however, at this stage 
of the analysis they lead to an additional uncertainty because the corresponding values of 
parameters do not well determined p4 |. Nevertheless, note here that in the framework of the 



analytic approach the phenomenological role of "non-perturbative" power corrections, which 
are controlled by the operator product expansion, is changed, and it will, in the future, be 
interesting to consider this problem further. 
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